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Abstract 


Frequency-domain methods are used to extract the open-^loop dynamics of the XV-15 tilt- 
rotor aircraft from flight test data for the cruise condition (V = 170 knots) The frequency 
responses are numerically fitted with transfer-function forms to identify equivalent modal 
characteristics The associated handling quality parameters meet or exceed Level II, Cate- 
gory A, requirements for fixed-wing military aircraft Step response matching is used to 
verify the time-domain fidelity of the transfer-function models for the cruise and hovet flight 
conditions. The transient responses of the model and aircraft ate in close agreement in all 
cases, except for the normal acceleration response to elevator deflection in cruise This dis- 
crepancy is probably due to the unmodeled rotor rpm dynamics The utility of the frequency- 
domain approach for dynamics identification and analysis is clearly demonstrated 


1. Introduction 



The identification of XV-15 tilt-rotor dynamics (Fig 1) from flight test data is an 
extensive ongoing effort to support the development of the next generation of tilt rotors - the 
joint services V/STOL aircraft (JVX) . 

The key concerns of the effort are the 
documentation of open-loop XV-15 
dynamics, and the validation of 
generic tilt-rotor models (Refs [1] 
and [2]) A frequency domain-based 
identification approach was developed 
and successfully applied for the hover 
flight condition (Ref [3]) Transfer 
function models describing the open- 
loop response characteristics of the 
XV-15 aircraft were extracted and com- 
pared with the simulation characteris- 
tics. Reference [3] presents a 
detailed description of the frequency- 
domain methodology and the results for 
the hover flight condition. 


As in all identification 
research, a key concern of this effort 
is the fidelity of the extracted 
models for input forms other than 
those used in the identification pro- 
cess The flight testing technique of 
Ref [3] uses a pilot-generated swept 
sine wave input to excite the vehicle 
dynamics. This yields an excellent 
identification of frequency responses 
and transfer functions for sinusoidal- 
like inputs. Such information is 
useful for frequency domain-based 
handling quality specifications such 
as the Mil Handbook for military air- 
craft (Ref. [4]) However, for time- 
domain specifications such as 
Mil-H-8501 (Ref [5]) and Mil-F-83300 
(Ref. [6]), criteria are largely based 
on responses to Step inputs, so 
extracted models must accurately 
reflect these characteristics as well 


Figure 1 The XV-15 Tilt-Rotor Aircraft (a) Cruise 
configuration, (b) Hover configuration, 
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This paper reviews the frequency-doinain based methodology and discusses the identifica- 
tion of transfer function models for the cruise flight condition (170 knots). This is a good 
limiting case for comparison with the previous hover results (Ref. [3]). Time-domain matching 
is presented to verify the step response characteristics of the extracted models for the hover 
and cruise flight conditions. These results show the utility of relatively simple transfer 
function models for handling qualities and control system applications. 


2. Review of Identification Methodology 

This section reviews the frequency-domain identification and time-domain verification 
techniques. The details of the identification approach are extensively discussed in Ref [3] 
and are only outlined in this paper. 

Frequency-domain identification is based on the spectral analysis of input and output 
time histories using Fast Fourier Transform techniques, This analysis produces describing 
functions which are (complex- valued) linear descriptions of the input-to-output processes The 
identification results are presented in Bode plot format: magnitude and phase versus frequency 

System bandwidth and effective time delay, important metrics in current handling qualities 
specifications, can be read directly from these plots Tabulated frequency response results 
are fitted with analytical transfer function forms to extract modal characteristics for handling 
qualities specifications given in terms of lower order (equivalent) system models Also, the 
transfer-function models can be driven with step inputs to extract familiar time domain metrics 
such as rise time, overshoot, and settling time 

The swept sine wave (frequency-sweep) input is a good excitation for the frequency-domain 
identification approach This excitation results in bounded and reasonable excursions of the 
aircraft, suitably exciting the important rigid body modes over the entire frequency range of 
interest (0, 2-6.0 rad/sec) The input is generated by the pilot in one axis at a time, with 
minimal regulation of the remaining axes. Starting from a trim condition, the pilot first exe- 
cutes two 20 sec period inputs to ensure good low-frequency identification Then the frequency 
of the inputs is slowly increased up to a maximum of about 6 rad /sec, yielding a total run 
length of about 90 sec Repeat runs are executed to allow concatenation of the 90 sec records 
in the analysis; this technique minimizes the effects of random noise 


When the bare airframe characteristics of the vehicle are highly unstable, as in the 
XV-15 in hover, execution of a 90 sec frequency sweep on the open-loop vehicle is clearly not 
practical. Subject to some important qualifications, the bare airframe transfer functions can 
be extracted from the closed-loop flight data Referring to Fig 2, identification of the open- 

loop transfer function (q/6g) requires 
(BARE' measurement of the output (q) and 

(PILOT) AIRFRAME) (OUTPUT) surface deflection (6^) signals This 

(INPUT) I total surface deflection is made up of 

components from the pilot (^lon^ 
stability and control augmentation sys- 
tem (SCAS) feedback (6f ) A key 
e 

requirement is that the total surface 
deflection must contain a significant 
component from the pilot which is uncor- 
related with the output The use of the 
"programmed" swept sine wave input 
avoids the otherwise troublesome result 
of identifying the inverse feedback 
transfer function (-1/Gp ) from passive 



tracking tasks 


(SCAS) 


Figure 2. 
model . 


Single deg ree-of -freedom pitch response 


When the vehicle dynamics are 
highly coupled, multiple control deflec- 
tions will occur in all degrees of free- 
dom. If the surface deflections are not 

fully QOVTelated with each other, the desired single-input transfer function can be isolated 
using the multi-input multi-output methodology described in Ref [3] This restriction was 
satisfied for the hover flight condition since the cross-control inputs supplied by the pilot 
were largely uncorrelated with the primary frequency sweep inputs, Therefore, it was possible 
to isolate the yaw response to ailerons- an important source of coupling for the hover flight 
condition. 


For the cruise flight condition, open- loop longitudinal transfer functions were also 
identified from the closed- loop tests since the elevator inputs supplied by the pilot were uncor- 
related with the outputs and the other surface deflections In the lateral axis, the large 
degree of inherent closed-loop stability allowed the pilot to execute frequency sweeps of the 
aileron without applying rudder corrections. However, the stability augmentation system supplies 
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rudder inputs which are highly correlated to the aileron inputs and roll rate outputs There- 
fore, extraction of the single- input, single-output open-loop transfer functions from closed- 
loop flight data was not possible for this flight condition; the lateral -directional tests were 
repeated with the SCAS disengaged, When the bare airframe dynamics are too highly coupled or 
unstable to make this practical, the open-loop dynamics can usually be extracted from the 
closed-loop test data, as was successfully accomplished for the XV-15 in hover (Ref [3]) 

Once the flight data have been digitized and preprocessed, the input, output, and cross 
spectra (G^xC^)* Gyy(f) , and GxyCf)* respectively) are calculated using modern Chirp 
z-transform methods (Ref [3]) Specific transfer functions G(f) are obtained from the ratio 
of the appropriate cross and input auto spectra: 


G(f) 


G (f) 
xy 

G (f) 

XX 


and are presented in Bode format The coherence function (Yxy) defined as 


""xy “ Tg 




XX* ‘ yy' 


( 1 ) 


( 2 ) 


is a good indication of the Input-to-output linearity This frequency-dependent parameter may 
be interpreted as that fraction of the output spectrum which can be accounted for by linear 
relation with the input spectrum, When the process under investigation is perfectly linear and 
the spectral estimates are noise free, the coherence function will be unity for all frequencies 
in the excited input spectrum range A value of the coherence function less than unity will 
result from nonlinearities in the system, input /output noise, or cross-coupled control inputs, 
The magnitude and phase responses are then fitted with analytical transfer function models to 
obtain closed form descriptions of the input-to-output processes In order to obtain a unique 
fit of the frequency responses, certain physical restraints on commonality of the transfer 
function denominator factors are imposed 


Current military handling quality specifications (Ref [A]) limit the allowable mismatch 
between the transfer function fits and the frequency response data, However, even small decibel 
deviations in the frequency-domain can produce surprisingly large discrepancies in time-domain 
correlations Time-domain verification first requires converting the transfer functions to 
canonical form: 


where 
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c = [t>0 ••• bj, 0 ••• o] 


and Sq, a^, • , are the denominator (ascending) coefficients and bg, bj^, . , b^ 

are the numerator (ascending) coefficients of the transfer function model being evaluated 
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This model is driven with step-input data from the flight tapes, and the resulting 
responses are then compared to the aircraft data for time-domain fidelity. This approach 
allows the verification of transfer functions associated with the stable degrees of freedom, 
since spurious model inputs (due to inexact initial condition matching and turbulence regula-^ 
tion) do not cause a divergence of the transients with time However, for the unstable degrees 
of freedom, these spurious inputs cause a rapid divergence of the flight and model responses 
For such cases, conclusive quantitative statements on the step response fidelity of the trans- 
fer function models cannot be made, even though the qualitative form of the response is correct 
and small differences in the predicted level of instability are probably not important 


3 Dynamics Identification for the Cruise Flight Condition 

In this section, the characteristics for the cruise flight condition are considered: 

V = 170 knots (indicated), nacelle incidence = 0®, altitude = 8000 ft The longitudinal and 
lateral flight dynamics are fully decoupled in this condition The primary longitudinal bare- 
air frame transfer functions of interest are pitch rate and normal acceleration responses to 
elevator deflection, q/6g and a^/Se* respectively The important lateral-directional transfer 
functions are roll rate response to aileron, p/6g, and sideslip (at the c.g, ) response to 
rudder, ^^g/^r* 

The following discussion of longitudinal dynamics identification and analysis is pre- 
sented in detail to illustrate the frequency-domain methodology and to compare the longitudinal 
results for the hover flight condition presented in Ref [3] 

Longitudinal Dynamics 

Since the longitudinal and lateral dynamics are fully decoupled and elevator surface 
deflections are measured directly, the flight testing for the longitudinal axis was conducted 
with the SCAS engaged This was not considered essential, however, since the dynamics for this 
flight condition are very stable 

The longitudinal stick displacements for two concatenated frequency sweeps are shown in 
Fig. 3 As previously discussed, each run lasts roughly 90 sec and is initiated with two 20~see 

cycles to ensure good low-frequency iden- 
tification The sinusoidal stick deflec- 
tion is very regular with a nearly con- 
stant amplitude of roughly ±5% The 
elevator signal is also very regular with 
a nearly constant amplitude of 2 deg 
These data are very typical of the 
transition-flight and forward-flight 
condition results. The relative stabil- 
ity of the longitudinal and lateral axes 
made execution of the frequency sweeps a 
very rapid and acceptable flight-test 
procedure The input auto spectrum for 
the concatenated elevator time histories 
is shown in Fig 4 This spectrum is 
fairly flat at low frequencies, dropping 
off slightly for the higher frequency 
inputs For frequencies outside the 
plotted range, the magnitude drop-off is 
much more pronounced 

The pitch-rate signal (conditioned 
with a 2 5 Hz low-pass filter for output 
presentation only) is shown in Fig, 5. 

The pitch-rate excursions are very regu- 
lar with a roughly constant peak-to-peak 
amplitude of 5 deg /sec, These oscilla- 
tions are eonsidered very acceptable to 
the research pilots, and are comparable 
with those encountered during the hover 
flight test The associated output auto- 
spectrum for the concatenated pitch-rate 
time histories are shown in Fig. 6. The 
output spectrum follows the input spectrum 
at low input frequencies and exhibit the 
familiar 20 dB/decade roll-off at higher 



Figure 3, Two longitudinal stick frequency sweeps 
cruise. 



Figure 4 Elevator surface input autospectrum 
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frequencies (which is typical of rigid 
body motion) . The cross-spectrum 
between elevator surface and pitch 
rate is shown in Fig> 7 and demon- 
strates a similar characteristic. 

The elevator surface-to-pitch 
rate transfet function, q/S^, is shown 
in Figs, 8(a) and 8(b) The smooth 
spectral data over the majority of the 
frequency range are typical of the 
excellent results obtained from the 
Chirp z-transform algorithm. Due to 
the short period-mode excitation, 
there is a peak in the pitch-rate 
response at abou^ 2-0 rad/sec, with 
an associated phase lag of 45 deg At 
higher frequencies, the magnitude 
response rolls off at 20 dB /decade and 
the phase shift approaches -90 deg 
owing to the K/s pitch-rate charac- 
teristic and negligible flexibility/ 
servo-lag effects. The drop in the 
magnitude response and associated 
positive phase response for frequen- 
cies below 0,3 rad/sec are due to the 
phugo id dynamics , 



Figure 5 Pitch-rate response (q) during longitudinal 
frequency sweeps 



The coherence function shown in 
Fig 9 is strong over the frequency 
range of interest (0 .2-6 ,0 rad/sec) , 
as shown. For input frequencies above 
7 0 rad/sec, the coherence function 
becomes erratic and the transfer- 
function identification is less accu- 
rate For low-frequency inputs (less 
than 0 3 rad /sec) , the pitch-rate 
response decreases, even for the 
nearly constant input amplitude owing 
to the effect of the phugo id dynamics 
This results in a decrease of informa- 
tion transfer and an associated drop 
in the coherence function. This 
coherence -function roll-off is also 
attributable to atmosphere turbulence 
effects, which become more important 
at low frequency where the turbulence 
spectrum peaks (Ref [7]) 


Figure 6 Pitch-rate output autospectrum (0^^) 



Figure 7 Magnitude of cross-spectrum between ele- 
vator surface inputs and pitch-rate response (iGfi^ql^ 


The Bode plot of Figs 8(a) and 8(b) is the basic format for presenting classical 
frequency-domain results It provides useful information on dynamic characteristics and on 
handling qualities implications To the extent that the linear-describing function is an accu- 
rate representation of the input-to-output process, the Bode plot is an exact description of the 
vehicle dynamics and contains no implicit assumptions on model order or structure, This is the 
key advantage of frequency-domain identification over the more conventional time-domain 
approaches 


Recent military specifications for piloted-handling qualities (Ref [4]) are based on 
two key frequeney-domain parameters: bandwidth (wg^) and effective time delay (xp) For a par- 

ticular transfer function, the bandwidth is defined as that frequency where the phase margin is 
45 deg, or the gain margin is 6 dB, whichever is lowest, As discussed in Ref [4], the band- 
width is a measure of the speed of response; a high bandwidth reflects quick response and accu- 
rate tracking capability; a low bandwidth suggests sluggish response and pilot- Induced oscilla- 
tion tendencies. The effective time delay is a measure of the slope of the phase curve for 
frequencies near the bandwidth value A large effective time delay (e .g , greater than 
100 msec) results in a significant degradation of piloted handling qualities. 


For longitudinal control, a key transfer function is pitch attitude response to elevator, 
6/dg, which is easily derived from the pitch-rate transfer function as; 


Q(s) ^ q(s) 
6^(s) ~ s6^(s) 


( 4 ) 
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Figure 8 pitch-rate response to elevator (q/-6^) 
in cruise (a) Transfer function magnitude, 

(b) Transfer function phase 



Figure 9 Coherence function for pitch-rate 
response identification (Y|^q) 


FLIGHT DATA 

TRANSFER FUNCTION MODEL 



Figure 10. Vertical-acceleration response to 
elevator (a^/S^) in cruise, (a) Transfer-function 
magnitude, (b) Transfer-function phase. 


where s is the Laplace operator 
The Bode plot for pitch-attitude 
response may be obtained from the 
pitch-rate response by rotating the 
magnitude curve about the m « 1 rad/sec 
point by slope of 20 dB/decade. The 
phase curve is shifted by a constant 
value of -90 deg. Referring to Fig. 8, 
the critical bandwidth criterion is 
phase margin which is 45 deg at a fre- 
quency of « 25 rad/sec The 

effective time delay is negligible since 
the phase curve is nearly flat where it 
is near 180 deg. These bandwidth and 
time delay values are well within the 
Ref [4] Level II, Category A, specifi- 
cations (adequate handling qualities for 
high-precision tracking tasks) , and cer- 
tainly reflect satisfactory character- 
istics for SCAS failure conditions. 

The transfer function for normal 
aaaeleration response to elevator^ 
a^/^e* is shown in Figs. 10(a) and 10(b) 
The response is dominated by the classi- 
cal second-order short -period mode over 
the majority of the frequency range. 

The magnitude curve is flat at mid- 
frequency indicating a constant normal 
acceleration response to a step elevator 
input, with a roll-off in response for 
frequencies beyond the short period mode 
The fall -off in normal acceleration 
response for frequencies below 
0 3 rad/sec is probably due to the 
dominance of the phugoid dynamics The 
phase curve exhibits the classical second 
order response; i.e , zero phase lag at 
low frequency, 180 deg of phase lag at 
high frequency, and 90 deg of phase 
lag at the second order mode 
(o3 ^ 2 ,0 rad/ sec) This is consistent 
with the previous pitch-rate results u 
The coherence function for the normal 
acceleration response shown in Fig 11 
is strong over the frequency range of 
0 3 to 10.0 rad/ sec, with the fall-off 
at low frequency, again owing to the 
dominance of phugoid dynamics and tur- 
bulence effects As before, this sug- 
gests excellent identification over the 
entire frequency range of interest. 

Longitudinal Transfer Function Fitting 

Analytical transfer function 
forms are selected for each degree of 
freedom based on configuration and flight 
condition factors, In the hover flight 
condition, vehicle dynamics are dominated 
by the hovering cubic, and decoupled 
heave and yaw modes. In wing-borne 
flight, the conventional longitudinal 
and lateral quartic equations dominate. 
Thus, transfer-function models which are 
appropriate for the hover flight condi-- 
tion are not necessarily applicable to 
forward flight conditions. Obviously, if 
a model of high enough order is selected, 
the parameters can be adjusted to accom- 
modate each flight condition. However, 
such transfer function models no longer 
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retain the physical significance of the 
classical lower-order parameters. 

Also, higher-order transfer-function 
models tend to be strongly tuned to the 
specific inputs which are used in the 
identification procedure (e .g , a fre- 
quency sweep) and are often very poor 
predictors of other test inputs (e g , 
step inputs) Therefore, higher-order 
models are not desirable. In the 
approach taken in Ref [3], the 
minimum-order transfer-function models 
which can satisfactorily fit the fre- 
quency responses are used with the 
upper limit taken as the physical order 
of the system. 



Figure 11, Coherence function for vertical- 
acceleration response identification 


Examination of the longitudinal frequency responses of Figs. 8 (a) and 8 (b) shows that the 
longitudinal dynamics of this flight condition are dominated by the short period mode 
(o)sp = 2 rad/sec) Therefore, we adopt the classical pitch rate and normal acceleration 
responses to elevator*: 


q(s) ^ ^ 

[ Cgp 9 ^^Sp3 


(5) 


where q(s)/ 6 g(s) is the Laplace-transformed pitch rate response to elevator surface deflection, 
deg/sec/deg-elevator ; is the elevator pitch sensitivity ^ the first-order numera- 

tor inverse time constant ^ tgp, Wgp are the equivalent short-period mode damping and natural 
frequency, respectively; and tq is the effective time delay; and 


-Tg s 

V « Z 

az(s) ^ ^ 

5e(s) [Cgp. “spl 


( 6 ) 


where a 2 (s)/ 6 g(s) is the Laplace-transformed vertical acceleration response (positive downward) 
to elevator surface deflection, g/deg-elevator ; is the elevator vertical sensitivity; and 

the denominator parameters are identical to those of Eq (5) The effective time delay for the 
vertical acceleration response is 

We ignore the low-frequency phugoid dynamics since, as seen from Figs 8 and 10, these 
are important only for the very lowest frequency inputs 

The longitudinal transfer functions of Eqs (5) and ( 6 ) have the same denominator factors 
so simultaneous fitting of these responses, while imposing the restriction of commonality of 
denominator parameters, makes good physical sense Also, the analyses of Refs, [3} and [ 8 ], 
suggest that such simultaneous fitting techniques are needed to ensure unique and physically 
meaningful values for the numerator factors For the longitudinal transfer functions, this 
avoids the otherwise so-called ’’fixed or free” L^ problem (Ref [4]) 

The transfer-function parameters of Eqs (5) and ( 6 ) are iteratively varied in order to 
obtain a best least-squares fit between the equations and the frequency responses of Figs 8 
and 10 over the selected frequency range This procedure is completed using the computer pro- 

hONFIT developed by Hodgkinson et al (Ref [9]) For the present case, the selected range 
for simultaneous fitting of the pitch rate and normal acceleration responses is 0 ,3 to 7 rad/sec 
In this frequency range, the dynamics are clearly dominated by the short-period mode and the 
coherence is strong for both transfer functions Once the short-period damping and natural fre- 
quency (Cgpj i^sp) obtained for the simultaneous fit, the high-frequency gain for the normal 

acceleration response (Zg^) is varied holding the damping and frequency constant. This opti- 
mizes the fit over the frequency range 03 to 10.0 rad/sec since the coherence function for 

this measurement remains strong out to higher frequencies The transfer function parameters for 
the pitch rate and normal acceleration responses are finally obtained as: 


^Shorthand notation: [C, w] implies s^ + 2^o)S + o)^, t = damping ratio, m = undamped 

natural frequency (rad/ sec) ; and (1/T) implies s 4- (1/T) , rad/sec 
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~ deg/sec^/deg-elevator'^ 

= -32 ,144 deg/sec^/in .~6^Qjg I 


I/T 9 = 1.035 rad/sec 
^ 2 



1 

Csp = 0 554 


= 2 179 rad/sec 

sp 


Tq = 0.016 sec 


= 1 597 g/deg-elevator"^ 


= 6 644 g/in - 6 xx)N 


Csp = 0 554 i 

> 

Wgp = 2 179 rad/sec 


T- =0 .018 Sec 
az J 



( 7 ) 


( 8 ) 


A comparison of these lower-order models with the flight test results is presented in Figs 8 
and 10 The pitch-response model matches the data very well over the selected fitting range 
(0,3“10,0 rad/sec for q/ 6 ^), which shows that the short period approximation adequately repre- 
sents the high- and mid-frequency dynamics For low- input frequencies (oj < 0,3 rad/sec), the 
phugoid mode causes a drop in the pitch response, a characteristic which is not ’’captured" by 
the short period model The match between the normal acceleration transfer function and the 
flight data is not nearly as good in the fitting range, and shows a noticeable discrepancy in 
the magnitude response for frequencies of greater than 2 rad/sec Possible sources of the mid- 
and high-frequency discrepancy are the rotor rpm dynamics, i^ich are known to be important for 
the forward flight conditions (Ref [10]} and rotor inflow dynamics, Neither of these sources 
of added dynamics is directly accounted for by the classical short-period approximation, except 
to the extent that the short period parameters are adjusted in order to best fit the flight 
data. As in the pitch response, the low frequency mismatch is due to the omission of the phu- 
goid mode from the transfer- function model 

Guidelines are given in the Mil Handbook (Ref [4]) for the maximum allowable mismatch in 
lower-order system modeling The discrepancies shown in the normal acceleration responses 
(Figs 10(a) and 10(b)) are well within these allowable guidelines, so the identified parameters 
should be useful indicators of piloted-handling qualities The capability to directly evaluate 
the effect of model order on response fitting is clearly seen in this example, again represent- 
ing a unique advantage of the frequency domain approach Further analyses could include the 
rotor rpm and phugoid degrees -of -freedom in the transfer- function models, and could evaluate 
their importance with respect to vehicle dynamics and handling qualities 

Lower order system parameters are useful in characterizing the vehicle dynamics for com- 
parison with the dynamics of other aircraft The identified short period damping and frequency 
(Csp “ 0.554 and ojsp = 2 18 rad/sec) are consistent with preliminary observations which are 
based on the raw Bode plot information. The relatively high degree of damping is a reflection 
of the small response peak of the pitch-rate transfer function The small effective time delays 
(in the pitch attitude and normal acceleration transfer fits) indicate that the high frequency 
flexibility and servo-lag effects are not important for this flight condition, as noted earlier. 
The Mil Handbook requirements for lower-order equivalent system pitch response are given in 
terms of the parameters WspT^^, Cgp, and T 9 The values of these parameters given in Eqs. (7) 
and ( 8 ) are well within the Level I, Category A, handling qualities requirements, confirming 
that the longitudinal characteristics for SCAS off conditions are adequate. 

Time-domain verification of the longitudinal transfer function models will be addressed 
after lateral response models are presented 

Lateral Dynamics 

The open-loop transfer function for voll-vate response to ailerons p/5a» is shown in 
Figs, 12(a) and 12(b) This is a classical first-order system in which aileron inputs produce 
a constant roll acceleration at high frequency, and a constant roll rate at low frequency.. The 
corner frequency is about 1 rad/sec, with an associated phase lag of roughly -45 deg, as 
expected, The dominant time constant is thus roughly 1 sec, implying about 2 to 3 sec to reach 
a steady-state roll rate. The coherence function shown in Fig 13 is very strong over the fre- 
quency range of 0 1 to 90 rad/sec, indicating an excellent identification of the roll-rate 
dynamics. 
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Specifications for the lateral 
handling qualities of military air- 
craft are g^.ven in the Mil Handbook 
(Ref [4]), based on a lower-order 
equivalent-system model for the roll- 
attitude dynamics (4>/d^) When roll/ 
yaw coupling is minor, the specifica- 
tion on the equivalent system roll 
mode time constant (1^) can be inter- 
preted as the reciprocal of the roll 
response bandwidth . Referring 

again to Fig. 12(b), the bandwidth for 
roll att'ttude response to ailerons is 
^BW<^ = 0 9 rad/sec, yielding a domi- 
nant time constant of about 1 .1 sec 
This satisfies the Ref [4J specifica- 
tion (Tj- <1.4 siec) for Level I, Cate- 
gory A, handling qualities 

Aerodynamic sideslip measure- 
ments are obtained from a sideslip 
indicator located roughly 18 ft ahead 
of the aircraft c g The sideslip at 
the c .g. , 3cg> is calculated by cor- 
recting the measured signal for posi- 
tion error based on yaw rate and air- 
speed No corrections are made for 
sensor dynamics because these are not 
fhlt to be significant within the 
ident if icat ion bandwidth 

The transfer function for side- 
slip response (at the c ,g, ) to rudder 
inputs Bcg/'^r* i® shown in Figs 14(a) 
and 14(b) The response is character- 
ized by a lightly damped second-order 
mode with a frequency of about 
1 6 rad/sec High-frequency rudder 
inputs yield a constant sideslip accel- 
eration, while low-frequency rudder 
inputs yield a constant sideslip angle 
corresponding t j a steady-state yaw 
rate The coherence function for this 
case, shown in Fig 15 is strong over 
the frequency range 0,1 to 5-0 rad/sec, 
falling off sharply in the high fre- 
quency range 

The rapid and relatively low- 
frequency decline of the coherence 
function of Fig 15 is unlike the pre- 
vious transfer functions considered 
for the cruise flight condition The 
reason for this difference is that 
aerodynamic sideslip is a state vari- 
able which is one derivative lower 
than the other rate variables consid- 
ered previously The response of 
parameters such as roll rate and 
pitch rate falls off at high fre- 
quency with a slope of 20 dB/ decade, 
However, the aerodynamic sideslip 
response as seen in Fig. 14(a) falls 
off at the much faster rate of 
40 dB/decade Although the input 
autospectrum is roughly constant 
over the entire plotted frequency 
range, this rapidly falling response 
causes a marked reduction of the 
output autospectrum and cross -t 
spectrum and an associated decline 




Figure 12 Roll-rate response to aileron (p/6a) in 
cruise (a) Transfer-function magnitude 
(b) Transfer-function phase 



Figure 13 Coherence function for roll-rate response 
identification (vl^p) 



Figure 14 Sideslip response to rudder (3cg/5j,) in 
cruise (a) Transfer-function magnitude, 

(b) Transfer-function phase 
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in the coherence function, This empha- 
sizes the advantages of using the 
higher state derivatives in identify- 
ing the high-frequency dynamics and 
suggests the possible advantages of 
using attitude signals to identify the 
low-frequency dynamics However, since 
the yaw mode of interest has a natural 
frequency of roughly 1 .6 rad/ sec, the 
falling coherence for frequencies 
greater than 5 rad/sec is not a severe 
limitation. 

As in the roll axis, handling 
qualities specifications for the yaw 
axis response are given in terms of 
equivalent system models These are 
derived in the following section. 


Lateral Transfer Function Fitting 

Equivalent system fitting using decoupled and coupled lateral/directional models has 
been considered in detail by Bischoff (Ref [8]) With the decoupled model, the responses are 
fit independently, using a first-order roll rate and a second-order sideslip transfer function 
The coupled model approach is based on simultaneous fitting of the roll and sideslip responses 
to obtain the following fourth-order transfer functions- 


P(s) 

«a(s) ■ (l/Ts)(l/TpUj. 


(9) 


where p(s)/6a(s) is the Laplace-transformed roll rate response to aileron surface deflection, 
deg/sec/deg-aileron; L§^ is the aileron roll sensitivity; are the second-order numera- 

tor damping and natural frequency, respectively; 1/Tg, 1/Tj- are the equivalent spiral and roll 
subsidence modes; ^d* equivalent Dutch roll mode damping and natural frequency, 

respectively; and Xx is the effective time delay 


5 ^( 3 ) 


YSr U/Tgp (1/Tg3)e""6® 

(1/T^)(l/Tr)[;d, Wj] 


( 10 ) 


where 3cg(s)/5r(s) is the Laplace-transformed sideslip response to rudder surface deflection, 
deg/deg-elevator ; is the rudder sideslip sensitivity; l/T^^^, l/T^g’ ^^*^33 first- 

order numerator inverse time-constants; and the denominator parameters are identical to those 
of Eq. (9) The effective time delay for the sideslip response is The simultaneous fit- 

ting approach is consistent with that used in the longitudinal axis and allows the identifica- 
tion of the tAJrfj/wd "coupling effect" in the roll response (Eq (9) , which is important in the 
handling qualities assessment 

The coherence function results of Figs 13 and 15 show a satisfactory identification of 
both roll and sideslip responses in the frequency range of 0,1 to 5,0 rad/sec, Based on this 
fitting range, the computer program LATFIT (Ref [11]) was used to obtain the parameters of 

Eqs (9) and (10), simultaneously Since the coherence for the roll response is satisfactory 

over the entire range of 0 1 to 10 ,0 rad/sec, this degree of freedom was refit alone, holding 
the denominator factors constant at the values obtained from the simultaneous solution This 

procedure is analogous to that used in the longitudinal case and optimizes the values of the 

high-frequency gain and numerator parameters in the roll transfer function. The final results 
for the parameters of Eqs. (9) and (10) are: 

Lx *= 4 486 deg/sec^/deg-aileron ^ 


- 17,630 deg/sec^/in -6' 
= 0 313 

0 ). - 1 ,887 rad/sec 


LAT 


1/Tg = 0 063 rad/sec 
1/T^ = 1.090 rad/sec 


- 0 .248 


( 11 ) 
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( 11 ) 
(cont ) 


< 0 ^ « 1,581 rad/sec ' 
T± * 0.045 sec 


051 deg/sec/deg-rudder' 
« -0.408 deg/sec/in. 

1/To * 0.086 rad/sec 

1/To - 0,818 rad/sec 

1/Tq =47 946 rad/sec 
P3 

1/T_ = 0 .063 rad/sec 


1/T^ = 1 .090 rad/sec 


= 0 .248 


u)^ - 1 .581 rad/sec 


"3 


0 026 sec 


( 12 ) 


The lateral /directional transfer function fits are plotted as dotted lines for compari- 
son with the flight data of Figs 12 and 14 The matching of the roll-response magnitude 
(Fig 12) is excellent over the entire frequency range Only a slight anomaly in the phase 
matching is apparent for low frequency inputs The agreement between the sideslip fit and the 
flight data is also excellent (Fig 14) 


As a result of the simultaneous fitting of the roll and yaw responses, the transfer 
function parameters of Eqs. (9) and (10) have clearly retained their physical significance. 

The inverse roll mode time constant, 1/T^ « 1 .090 rad/sec is roughly equal to the roll response 
bandwidth as expected; both are well within the Level I roll response requirements, as mentioned 
earlier. The equivalent Dutch roll mode is lightly damped with a natural frequency roughly 
corresponding to the peak in the sideslip Bode magnitude plot (Fig 14(a)) The small effective 
time delays for the roll and sideslip responses reflect negligible high-frequency flexibility 
and servo-lag effects, and are consistent with the previous longitudinal results. The small 
time delay in the sideslip response also supports the omission of corrections for sideslip 
sensor dynamics. 


The ratio of the natural frequency of the numerator complex zero (wp to the denominator 
complex pole (w^j) is one measure of roll/yaw coupling in response to aileron inputs When 
“ ^d’ numerator and denominator quadratic factors roughly cancel, and the resulting 
decoupled roll response is characterized entirely by the roll-mode time constant This case 
leads to the best handling qualities for a nominal value of the roll-mode time constant. As 
the roll /yaw coupling increases, the numerator and denominator quadratic factors of Eq (9) no 
longer cancel and an undesirable oscillatory component of roll rate is generated. Referring to 
Eq (11), the near unity value = 1 19 suggests no such concern for roll /yaw coupling 

A more direct measure of the roll/yaw coupling is obtained from evaluating the ratio of 
roll attitude to sideslip evaluated at the Dutch roll frequency l<|>/3|(| Analyses show 
(Ref . [4]) that there is a direct correlation between the |<})/3id ^sitio and the ratio of oscil- 
latory roll rate to steady roll rate for step aileron inputs A small value of |4>/3|d (e g,, 
less than 1 .4) suggests a small oscillatory roll-rate component and desirable handling quali- 
ties. A large value causes roll response overshoot and related poor tracking characteristics 
A cross correlation of roll rate and sideslip responses yields l<l>/3lxi ® 1 27, which implies no 
such tracking deficiencies. The preceding analyses have shown that the roll response is largely 
decoupled and first order in nature. Therefore, the simple, first-order roll rate model could 
be adopted for future studies 


4 . Time Domain Verification of Transfer Function Models 


The transfer function is the minimum realization description of a linear input-output 
process. The transfer-function models developed above c^n be used to generate all other fre- 
quency and time domain information, Therefore, a set of verified transfer function models is a 
useful output format for the linear identification procedure. However, the adequacy of these 
linearized transfer-f unction models for predicting the time-domain dynamics of nonlinear systems 
is at question. For such nonlinear systems, the transfer function is actually a describing 
function which is strictly yalid only for the input amplitudes which were used in the flight 
test experiment. Also, since identification procedures tend to be tuned to the type of test 
inputs used (e.g. , a sine-sweep), their accuracy in predicting responses to other classes of 
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inputs (e,g., steps, doublets) is often uncertain. Lastly, although the mismatches observed in 
the previous section between the transfer-function fits and the frequency-response flight data 
never exceeded the allowable handling quality specifications, the importance of these discrep- 
ancies in predicting time response behavior is still at issue. These questions are addressed 
in the following section. 

The basic verification approach is to compare the step response of the aircraft with the 
response of the transfer-function model driven by the same recorded control inputs. The step 
responses were completed for all flight conditions with the stabilization system disengaged in 
order to show that the extracted transfer-function models are valid for the bare-airframe con- 
figuration, The transfer functions for the hover flight condition and those extracted in the 
present analysis for the cruise flight condition are summarized in Table 1. For illustrative 
consistency, all transfer functions are referenced to the three surface deflections: elevator, 

aileron, and rudder The verification results for the hover flight condition are presented 
first 


TABLE 1.- SUMMARY OF IDENTIFIED TRANSFER-FUNCTION MODELS^ FOR THE 
HOVER AND CRUISE FLIGHT CONDITIONS 


Hover (V - 0) 


Cruise (V =170 knots, indicated) 


Longitudinal dynamics 


-0.108 e"° 
(0,115) 


; g/in.- power lever 


1.597 e~° 


[0.554, 2 179] 


; g/deg- elevator 


_ -2.6Q7s(-0.075)(0.686)e"° 

(0 115) (1,59 7) [-0,978, 0 393] * 

deg/sec/deg ■ elevator 


± = ‘‘7. 727(1.035)e^° ° 
6^ [0 ,554, 2 179] 


deg/sec/deg - elevator 


Lateral dynamics 


p „ 3.959s(-0.092)(0.570)e"'' p 

6, ~ (0 108) (1 429) [-0 581, 0.427] ’ 6 

3 3 

deg/sec/deg- aileron 


r 


6 


r 


0.708 e~° 
(0 108) 


deg/sec/deg - rudder 


^_C8 

«r 


r „ 0.398(-0. 310) [0.762, 0.524]e~'’ 

" (0. 108) (1 .429) [-0 581, 0 427] 

deg/sec/deg- aileron 


4.486s[0 313, 1.887]e““ 

(0.063)(1,090)[0 248, 1.581] ’ 

deg/sec/deg’^ aileron 

-0.051(0.086) (0.818) (47.946)0-°*°^^^ 

(0. 063) (1 090) [0 248, 1,, 581] * 

deg/deg - rudder 


‘^Shorthand notation: [^, w] implies s^ + 2Co)S + C = damping ratio, o) = undamped natural 

frequency (rad/sec) ; and (1/T) implies s + (1/T) , rad/sec 


Hover Flight Condition 

A co-plot of the transfer-function model for yaw rate response to rudder r/6- (Table 1) 
with the frequency response data of Ref [3] is presented in Fig 16 The transfer-function 
model matches the flight data very well at high frequency in both magnitude and phase but not as 
well at low frequency Even so, the fit meets the mismatch criteria of Ref [4] for handling 
qualities assessment These criteria emphasize the piloting crossover range of 1—3 rad/sec, in 
which the fit is excellent for this case 

When the transfer-function model is driven with one of the swept stne’^wcwe tnputs used in 
the identification procedure, the time-domain ramifications of model mismtch are clearly 
exposed* The transfer-function model and flight-data responses for a typical frequency sweep 
are compared in Fig 17 The input amplitude shown in Fig. 17(a) is roughly ±20% of the maximum 
pedal travel producing a yaw response amplitude of roughly ±20 deg/sec, as seen in Fig. 17(b). 
The flat spot in the yaw rate signal at 50 sec is due to yaw-rate sensor limiting. In general, 
the match between the transfer-function response and the flight data is excellent. The phasing 
and form of the model response matches the flight data very well, with the only noticeable 
anomaly being an underestimation of the peak amplitudes for low-frequency inputs. This is a 
result of the fit mismatch at low frequency (Fig 16) Thus, the transfer function model accu- 
rately captures the essence of the yaw-rate response over the frequency range of interest, and 
can be considered adequate for sine -wave type inputs 
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The flight data for yaw rate 
response to a st^p pedat 'input is shown 
in Figs. 18(a) and 18(b) The control 
deflection represents roughly 40% of 
the maximum control power, a fairly 
large amplitude input. The maximum 
yaw-rate response is about 35 deg/sec, 
which is 75% greater than the peak 
values obtained in the frequency 
sweeps The transfer-function model 
response co-plotted in the dashed 
curve of Fig. 18(b) compares extremely 
well with the flight data over the 
entire run- even for this fairly large 
input amplitude case Analyses on the 
vertical response to collective 
(a^/^c) showed the same excellent 
response correlation. This is 
expected since the response dynamics 
for the heave and yaw axes are nearly 
the same (Table 1) 

The step-response run lengths of 
roughly 10 sec emphasize the high- and 
mid-frequency ranges, allowing only 
enough time for about one time constant 
(Ty =9.3 sec) Thus, the low- 
frequency mismatches never really have 
enough time to build up This is 
typical of step-response data and 
illustrates the advantages of a sym- 
metrical sine -sweep input Clearly, 
however, the response in the first few 
seconds is the key concern of the 
pilot This is reflected in the mis- 
match criteria being most stringent in 
the mid-frequency range. 

So far, the time response 
results show that the simple transfer- 
function models match the step-response 
characteristics very well Also, the 
use of a very low-order description has 
avoided any problems in over- tuned 
modeling. In general, however, good 
time history matching is highly depen- 
dent on the level of spurious inputs 
to the model due to; (1) improper 
setting of initial conditions, 

(2) control inputs required to sup- 
press turbulence, and (3) unmodeled 
control coupling When the model 
dynamics are stable^ the effects of 
these spurious inputs die out with 
time, and response matching presents 
no particular problems Such cases 
are the easiest for time-domain iden- 
tification as well 




Figure 16 Yaw-rate response to rudder (r/6j.) in 
hover (a) Transfer-function magnitude, 

(b) Transfer-function phase 


— FLIGHT DATA 

-- -- TRANSFER FUNCTION MODEL 




Figure 17 Comparison of aircraft and transfer- 
function model responses to a pedal frequency-sweep 
in hover (a) Rudder surface deflection (b) Yaw- 
rate response. 


When the dynamics of the transfer-function response are unstable ^ the transients due to 
these spurious inputs grow with time Therefore, small errors in the choice of initial condi- 
tions, or small inputs required to suppress the response to turbulence, can create gross differ- 
ences between the transfer- function response and the flight data. Also, for unsta.ble vehicles, 
acquisition and maintenance of a steady trim condition is often not possible, further complicat- 
ing the proper initialization of the transfer function calculations. As the level of instabil- 
ity of the transfer-function dynamics increases, so do the problems of time-history matching. 

Not surprisingly, these are also the worst cases for time-domain identification. 

The roll-rate-to-aileron transfer function (Table 1) is dominated by an unstable roll 
oscillation having a time-to-double of about 2.8 sec This reflects the speed at \»hich the 
errors and inexact control inputs and initial conditions will propagate. A good step response 
in the roll axis was achieved starting from a fairly steady initial trim condition as shown in 
Figs. l9(a) and 19(b) The assumed initial condition is determined from the average of the 


J 
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slope, rise time, and percent overshoot 
are all well predicted; therefore, the 
transf er~f unction model is a useful 
characterization of the acceleration 
response This comparison further indi- 
cates the need to investigate sources 
of unmodeled dynamics (e.g., engine 
response) in the vertical axis. 

The aircraft roll rate response 
to a step aileron input is shown in 
Figs. 21(a) and 21(b) The transfer- 
function model shown in the dashed 
line accurately predicts that the 
response will be predominantly first 
order, with no overshoot or oscilla- 
tory tendency This corroborates the 
same conclusion made earlier on the 
basis of the small and 

ratios The sideslip response to a 
step rudder input is shown in 
Figs. 22(a) and 22(b) with the appro- 
priate corrections for the nose-boom 
position The transfer-function model 
shown in the dashed curve matches the 
flight data fairly well for the first 
10 sec. Thereafter, the responses 
diverge, probably due to a gust 
encounter or some other aerodynamic 
interference effect at the sideslip 
vane during the recovery phase Once 
again, the character of the step 
response is well predicted, with an 
accurate modeling of the rise time, 
overshoot, and steady state sideslip 

This verification study shows 
the utility of fairly simple transfer- 
function models in predicting the 
transient responses to relatively 
large and varied inputs High- 
frequency and mid-frequency spectral 
response matching is clearly essential 
for accurate step response prediction, 
as it is for handling qualities 
assessments, Time response matching 
for unstable transfer functions is 
very difficult because of the inabil- 
ity to acquire and maintain a steady 
initial trim, and the tendency toward 
rapidly growing transients due to 
small unmodeled excitations 



Figure 20. Comparison of aircraft and transfer- 
function model responses to a step elevator input in 
cruise (a) Elevator surface deflection, (b) Pitch- 
rate response (c) Normal (positive downward) accel- 
eration response 

— FLIGHT DATA 

TRAIMSFER FUNCTION MODEL 




Figure 21 Comparison of aircraft and transfer- 
function model responses to a step aileron input in 
cruise (a) Aileron surface deflection, (b) Roll- 
rate response, 
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Figure 22. Comparison of aircraft and transfer-function model responses to a step rudder input 
in cruise. (a) Rudder surface deflection (b) Sideslip response (at the c.go). 


5. Conclusions 


Some specific conclusions from this study of the open-loop dynamics of the XV-15 tilt- 
rotor aircraft are; 

1) The response characteristics for the cruise flight condition are very stable and 
decoupled r The handling quality parameters meet or exceed Level II, Category A, requirements 
for fixed-wing military aircraft 

2) Standard lower-order equivalent models adequately match the identified cruise dynam- 
ics in all degrees-of-freedom except for the normal acceleration responses to elevator deflec- 
tion The mismatch for this response is probably due to the unmodeled rotor rpm dynamics. 

3) The step responses of the identified transfer functions and the aircraft generally 
match very well for both the hover and cruise flight conditions. Time-domain verification of 
the unstable attitude dynamics in hover is very difficult because small errors at the beginning 
of the run cause the responses of the model and the aircraft to diverge rapidly. 

The frequency-domain approach has proven to be a relatively simple and accurate means 
for extracting the bare-airf rame dynamics of the XV-15 The utility of the derived lower order 
transfer-function models for handling qualities and control^-system studies has been shown in 
the frequency and time domains. 
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